function [Di Xi] = Interpolate(D, X)

% Dubbele tijdstippen verwijderen..
[B,I,J] = unique(X);
% I bevat locaties van de dubbelen.
Xtf = false(size(X));
Xtf(I) = true;
% Xtf
% find (~Xtf)
% X(~Xtf)
for i = X(~Xtf)
    Xtf(X==i) = false;
end
% Data cleaned from doubles..
Dc = D(Xtf);
Xc = X(Xtf);

% Stel hier de gewenste X resolutie mee in..
Xi = linspace(0,1,1000);
Di = interp1(Xc, Dc, Xi, 'lineair', 'extrap');

% Hoeft niet meer vanwege 'extrap'
% n = isnan(Di);
% Di = Di(~n);
% Xi = Xi(~n);

end